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Photons cannot be localized in a sharply defined region. The expectation value of their energy 
density and the photon number density can only be approximately localized, leaving an exponential 
tail. We show that one may sharply localize either electric or magnetic (but not both) footprints 
of photons, and only momentarily. In the course of time evolution this localization is immediately 
destroyed. However, the coherent states, like their classical counterparts, can be localized without 
any limitations. The main tool in our analysis is a set of space-dependent photon creation and 
annihilation operators defined without any reference to the mode decomposition. 

PACS numbers: 42.50.-p, 03.65.Ta, 14.70.Bh 



I. INTRODUCTION 

Physicists have pondered over the problem of pho- 
ton localization and the related problem of the photon 
wave function for almost 80 years now, beginning with 
the work of Landau and Peierls An extensive re- 
view of the photon localization problem was recently pre- 
sented by Keller . The problem of photon localization 
is closely related to the widely studied problem of the 
photon position operator. In this paper we introduce 
an operational definition of partial localization, based on 
the measurements of correlation functions for electric or 
magnetic fields. For a want of better names, we shall use 
the terms electric localization and magnetic localization 
even though this might erroneously suggest the presence 
of some electric or magnetic devices that confine the pho- 
tons. Since a sharp localization of photons according to 
our operational definition of localization is not possible, a 
photon position operator compatible with this definition 
does not exist. 

Earlier studies of the photon localization ernphasized 
the mathematical aspects (see, for example, [3, 0, [H])- 
In this paper we emphasize the physical properties of 
the electromagnetic field. In particular, we exhibit the 
role of the photon helicity and the symmetry between 
the electric and magnetic fields. We proceed in the foot- 
steps of Glauber 00, H, Q who was the first to recog- 
nize the significance of the space-dependent creation and 
annihilation operators. This approach was recently sum- 
marized and expanded in an extensive paper by Smith 
and Raymer [lO|. In our work we concentrate on the 
analysis of the photon localization in terms of the elec- 
tric and magnetic field operators. We treat these fields 
(after smearing over space-time regions) as bona fide ob- 
servables. We put emphasis on the field aspect that is 
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complementary to the particle aspect. It might have a 
weaker connection with experiments (usually based on 
photon counting as highlighted by Glauber) but it is more 
precise as was explained in detail by Bohr and Rosenfeld 
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The standard method of quantization of the free 
electromagnetic field, based on the decomposition into 
monochromatic modes, is not well suited for the discus- 
sion of localizability because the monochromatic mode 
functions are not localized. To overcome this problem we 
further developed an alternative method of quantization 
that does not require a mode decomposition. The essen- 
tial mathematical tools in our analysis are the Riemann- 
Silberstein (RS) vector and the helicity operator. This 
formulation has some merits of its own and it can also 
be used to study other general properties of photons (for 
example, the local aspects of entanglement) without re- 
strictions due to a specific choice of modes. 

Photons are quanta of the electromagnetic field — they 
are the carriers of the electromagnetic field hence every 
photon has its electric and magnetic side. The localiza- 
tion of photons studied in this paper relies on the elec- 
tric or magnetic manifestations of the photon's existence 
— on the footprints that a photon leaves. We say that 
the photons are electrically or magnetically localized in a 
region R if their electric or magnetic properties are con- 
fined to the localization region. The electric localization 
is complementary to the magnetic localization — elec- 
trically localized states are magnetically delocalized and 
vice versa. Since the electric and magnetic field vectors 
are canonically conjugate variables, this is akin to the 
complementarity of position and momentum in quantum 
mechanics. States of particles sharply localized in posi- 
tion are spread out in momentum and vice versa. 

The localization of photons discussed in this paper is 
not easily seen in the standard formalism since the mode 
functions that are introduced there correspond to well 
defined frequencies and one needs their superpositions to 
produce localized states. A much more convenient tool 
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is a set of space-dependent creation and annihilation op- 
erators discussed in Sec. IIIII These mode-independent 
operators enable us to introduce in Sec. IIVI the Fock 
bases which give a precise meaning to photon states. In 
Sec. I VII we introduce the notion of electrically or mag- 
netically localized photon states and in Sec. I VIII we show 
that their localization is an ephemeral effect — it is im- 
mediately destroyed by the time evolution. A weaker 
form of localization: diffuse localization is discussed in 
Sec. IVIIII and finally in Sec. IIXI we analyze localization 
in terms of coherent states. 

Throughout this paper we shall work in the Heisenberg 
picture. The observables will be represented by time- 
dependent operators while the state vectors will be con- 
structed at a reference time, taken as usual to be t = 0. 
Our definition of the creation and annihilation operators 
does not use the Hamiltonian or the time-dependence of 
the field operators. We do not split the field operators 
into its creation and annihilation parts according to the 
sign of the frequency. Instead, we rely solely on the equal- 
time canonical commutation relations. 



II. RS VECTOR AND THE PHOTON 
HELICITY OPERATOR 



lowing matrix: 



X = ns 



—iuz iuy 
in^ -in^ , , 
—iUy irix 



and is equivalent to the following cross product 



X 



iky. 



(4) 



(5) 



For plane waves, the positive (negative) helicity corre- 
sponds to the left-handed (right-handed) circular polar- 
ization. 

In position space, the cross product with the vector ik 
becomes the curl and the division by k is represented by 
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(27r)3 k " 27r2|r-r'|2' 



(6) 



Therefore, the action of the helicity operator x on ^-ny 
vector function V{r) becomes the following nonlocal op- 
eration: 



One hundred years ago Silbcrstein 

[HQ EKE 

[isj discovered that the following complex combination of 
the electromagnetic field vectors. 
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(1) 



is very useful in the analysis of solutions of the Maxwell 
equations. Indeed, the two pairs of Maxwell equations 
become one pair without any loss of information 



In the subspace of divergence-free vectors, the helicity 
operator is an involution, x" = 1. This property implies 
that the operators P±, 

P± = l^, P|=F±, P+P-=0, (8) 

are projectors on the positive and negative helicity sub- 
spaces. In Sec. IIIII we show that with the use of the 
helicity operator we can split the operator F{r,t) into 
its creation and annihilation parts. 



dtF{r,t) ^ -icV X F{r,t), V-F{r,t)^0. (2) 

One can always recover the D and B fields by taking the 
real and imaginary part, but treating F{r,t) as one ob- 
ject simplifies the analysis. In quantum electrodynamics 
the classical function F{r,t) is replaced by the operator 
F{r,t) in the Heisenberg picture satisfying the Maxwell 
equations. 

In particle physics the helicity operator x is defined as 
the projection s-n oi the spin operator s on the direction 
of momentum n = k/k. The helicity of photons, and 
other massless particles, has a Lorentz-invariant mean- 
ing. For spin-one particles the spin operator is a vector 
whose components are 3x3 matrices. The components 
of these matrices can be expressed in terms of the Levi- 
Civita symbol in the following way (in units of h): 



(3) 



Thus, the action of the helicity operator on any vector 
function in momentum space is represented by the fol- 



III. SPACE-DEPENDENT CREATION AND 
ANNIHILATION OPERATORS 

The components of the RS operator F{r,t) and its 
Hermitian conjugate satisfy the equal-time commutation 
relations (repeated indices always imply summation) 



F,{r,t),Fj{r\t)\ = hcea,di5^^\r - r'), (9) 



that follow directly from the canonical commutation re- 
lations [l^, [13] for the electromagnetic field operators 
D{r,t) and B[r,t), 



B,{r,t),Dj{r\t)] = -iheajdi5^^\r ~ r'). 



(10) 



By acting on the RS operator with the projection oper- 
ators 



(11) 
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we split F{r,t) into its annihilation, , and creation, 
JF"", parts. This identification follows from the commu- 
tation relations for these two parts of the RS operator 
(the remaining commutators vanish): 



F+ir,t),F+\r',t)'\ ^hcc±{r-r'), 
F- {r,t), Frt (r', t)] = -ftc c-^ (r - r') 

where the functions c,^ are 



(12a) 
(12b) 



The operators [ ([T5)) ] represent the annihilation and cre- 
ation parts of the field operators, 



D{r,t) = \l'^ [dt(r,t) + d(r,i)] , (18a) 



B{r,t) = \l^[b\r,t)^b{r,t)] 



(18b) 



The roles of these annihilation and creation operators 
will be fully elucidated in Sec. IIVI where we construct the 
Fock bases. 



(13) 



47r2 



r — r' 



These relations are obtained from the definition of the 
projector operators [([8])] and from the commutation re- 
lations [([9])]. The significance of this result is best seen 
in Fourier space, 



-zb/i ^ ^ijk^ kikj i 



(14) 



Both Hermitian matrices cjAk) and c-Ak) have the 
eigenvalues k and 0. The minus sign in the second com- 
mutation relation [((T^)] means that the operator F^ 
must be interpreted as a creation operator. In accor- 
dance with this interpretation, we introduce the following 
notation: 



(15a) 
(15b) 
(15c) 
(15d) 



d(r,t) = 


F+{r,t) 


+ F-''{r,t) 




/Wc 


d\r,t) = 


F-{r,t) 


+ F+^(r,i) 




/Wc 


b{r,t) = 


F+(r,t) 


-F-t(r,i) 




Jhc 


b\r,t) = 


F'{r,t) 


-F+'<{r,t) 




/he 



It follows from their definition that the operators d and 
b are not independent. 



b(r,i) = -ixd{r,t), 
b\r,t) = ixd\r,t)- 



(16a) 
(16b) 



All four vector operators d, d^, 6, and are divergence- 
free. Their commutation relations are 



d,(r,t),d](r',<)j = {d,d, - 5,, A) ^ (17a) 
J J zn \r — r \ 



2TT^\r — r 



6,(r,t),4(r',t)l = -leuA^^^K^-r')- 



-2^ (17b) 
(17c) 



IV. PHOTON NUMBER OPERATOR AND THE 
FOCK BASES 

We shall use the expression for the operator of the total 
number of photons, given by Zeldovich [23|. It has the 
following nonlocal form: 



N = 



1 



d'r I (fr' 



b{r,t)-D{r',t) , B{r,t)-B{r\t) 



e\r — r 



/|2 



/i|r — r 



l\2 



(19) 



We have omitted the time argument of N because the 
photon number operator is time- independent, as a result 
of Maxwell equations. With the use of Eqs. (fTS]) . we can 
explicitly carry out the normal ordering and obtain, 

J J 47r2|r-r'|2 

(20) 

One can prove, using Eqs. (|16p . that the contributions 
from the electric and magnetic parts to the total num- 
ber of photons are equal. Using formula (|A2p one can 
transform the number operator to the form, 

N= [d^r\ip^{r,t)-ip{r,t) + <f>Hr,t)-(i>{r,t)\ , (21) 



where the operators xp and cj) introduced by Cook [2S 
are connected with d and b through the equations 



V'(r,<) = 
4>{r,t) = 



1 



87r3/2 
1 



\r - r'|5/2 



87r3/2 J |^_^/|5/2 



d{r',t), 



b{r',t). 



(22a) 
(22b) 



The creation and annihilation operators satisfy the 
standard commutation relations with the photon num- 
ber operator; 



AT, d\rA) 




'N,b\r,t) 


^b^{r,t), 
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obtained with the use of the commutation relations ((T7)l . 
Eq. (|Aip . and the formula 



(24) 



As was to be expected, the creation operator increases 
and the annihilation operator decreases the number of 
photons by one. 

The photon creation and annihilation operators 
{r,t),b^r,t),d{r,t), and b{r,t) differ from their stan- 
dard counterparts: they are not labeled by quantum 
numbers (such as momentum and angular momentum) 
or by the mode number I enumerating some wave-packet 
modes Vi{r) as in 0, Instead, they are labeled by 
space points and they also have a vector character. We 
emphasize that in their definition we have neither used a 
mode decomposition nor the equations of motion. 

Having established the basic properties, we can pro- 
ceed further and construct the Fock basis. To this end 
we choose t — as our reference time for the construction 
of state vectors in the Heisenberg picture, and we drop 
the time parameter. We can construct the Fock basis 
with the use of either <i^(r) or 6^(r) operator. We be- 
gin with the operators dt(r). They furnish the space of 
states with the structure of an inverted pyramid of Fock 
states 



dl{r)dUr')dlir")4{r"')\0) n 



dlir)d]ir')dlir")\0) 

dlir)d]{r')\0) n^2 

4(r)|0) n^l 

|0) n = 



(25) 



These state vectors form a complete basis — every state 
can be described by the superposition of the basis vectors. 
Thus, an arbitrary n-photon state has the form 



dVi . . . / dV„(p*-*"(ri, . . . r„)4 (n) . . . <(r„)|0). 



(26) 



This basis is not orthogonal due to nonlocal commutation 
relations between the operators d{r) and <i^(r'). The sec- 
ond Fock basis can be constructed with the use of the op- 
erators 6^(r). These two bases are dual ^23*] with respect 
to each other (their vectors are mutually orthogonal) be- 
cause the commutation relations between b{r) and d^r') 
are local. In particular, the orthogonality condition for 
one-photon basis vectors has the form: 

(0|6,(r)4(r')|0) = -leuA^^^^r - r'). (27) 

The orthogonality for n-photon state vectors also follows 
from commutation relation (|17cp . 



V. ELECTRIC AND MAGNETIC PHOTON 
WAVE FUNCTIONS 

Every one-photon state can be obtained by acting on 
the vacuum state |0) with the creation operator con- 
structed from either d^ or b\ 



(frip{r)-d\r)^ d^r(fi{r)-b\r 



(28) 



where <p(r) and <p(r) are complex-valued vector wave 
functions. To guarantee the equality of both integrals, 
these two wave functions must be connected by the rela- 
tion 



(29) 



We shall assume throughout this paper that all vector 
wave functions are divergence-free. The normalization 
condition for both these functions has a nonlocal form, 



1 

2^ 



d^r 



d^r 



,[Vx<^*(r)]-[Vx v(r') 



= 1, 



(30) 



because of the nonlocal form of commutation relations 
(I17p . This condition guarantees the standard normaliza- 
tion of the commutator [a, a^] = 1. 

Even though both wave functions in Eqs. ((28|l can be 
used equally well to describe every photon state, one may 
be preferred over the other depending on the situation. It 
is natural to call ip{r) as the electric wave function and 
ip{r) as the magnetic wave function. Every particular 
wave function can be used in two different ways: it can 
be treated as an electric or magnetic wave function, lead- 
ing to two different photon states. These two states differ 
by an interchange of their electric and magnetic proper- 
ties (duality transformation). The existence of two pho- 
ton states described by the same function reflects perfect 
symmetry between electricity and magnetism of the free 
electromagnetic radiation field. 

The functions y(r) and <^(r) do not coincide with the 
photon wave functions studied before [l^, [13, [iM [111 but 
they are directly related. The photon wave functions 1F± 
for both signs of helicity were defined in as the matrix 
elements of the RS operator and its Hermitian conjugate 
taken between the vacuum and an arbitrary one-photon 
state vector. 



:P+(r,t) = (0|F(r,t)at|0), 
T-{r,t) = (0|Ft(r.,t)at|0). 



(31a) 
(31b) 



By adding and subtracting these two equations, with the 
use of Eqs. and we obtain at t = 



T+{r)+T-{r) ^ /^3^/^o|d(r)^(r')-bt(r')|0), 



he 



= — iV X ip{r) 



(32a) 



he 



= d^r'{0\b{r,0)ip{r')-d^{r')\0), 



V X ip{r). 



(32b) 
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These relations between both types of photon wave func- 
tions can be extended to all t because the wave func- 
tions obey the same evolution equations obtained from 
the Maxwell equations. 



VI. ELECTRIC AND MAGNETIC 
LOCALIZATIONS 

We adopt the following operational definition of the 
electric (magnetic) localizability. A photon state is elec- 
trically (magnetically) localized in some region R if all 
measurements of the electric (magnetic) fields (including 
their correlations) carried out outside of R do not reveal 
the existence of the photon. 

In order to construct localized photon states we choose 
a sharply localized (vanishing outside of R) vector func- 
tion ipii{r). Depending on whether we treat this function 
as an electric or a magnetic photon wave function, we ob- 
tain two different photon states created by the operators, 



namely. 



a)j, and a^. 



.ji- J (f'r'tpii{r)-d\r), 
4- Jd^riPR{r)-b\r). 
They obey the following commutation relations: 



r)(r,0),aj, 



B(r,0),a} 



r)(r,0),4 



-B(r,0),4 



^ V X ipRir), 



V X xipR{r). 



(33a) 
(33b) 

(34a) 
(34b) 
(34c) 
(34d) 



We shall show that these two operators create either mag- 
netically or electrically localized states. Indeed, when r 
lies outside of i?, the operator a|j commutes with the op- 
erator B(r,0). This property guarantees the magnetic 
localizability. To show this, let us consider the results 
of the measurements of the correlation functions of the 
magnetic field. They are described by the expectation 
values of any number of the operators of the magnetic 
field. 



(0|ai?,B,(ri,0)Bj(r2,0). 



.Bfc(^n>0)ajj|0). 



(35) 



When all points lie outside of i?, the operator a}j com- 
mutes with all field operators Bi{ri). Therefore, we may 
bring all the way to the left, and with the use of nor- 
malization condition (|30p we obtain the expectation value 
in the vacuum state, as if the photon was not present. 



\aRB,{ri,0)Bj{r2,0) . 
-- (0|B,(ri,0)4(r2,0) 



Sfc(r„,0)aj^|0) 
•Sfe(r„,0)|0). 



(36) 



According to Eq. (|34aP the correlation functions of the 
electric field involve the function x'ij}R{r). We prove in 
Appendix |B] that this function is not localized. There- 
fore, the expectation values of the products of electric 
field vectors in the magnetically localized photon state 
differ from the vacuum values outside of the localization 
region — the electric field is felt throughout all space. 
One may say, therefore, that in such a state the mag- 
netic properties are hidden inside of R while the electric 
properties are spread throughout the whole space. 
Creation operator (|33b[) can also be chosen in the form: 



(37) 



The one-photon state created by a}^ is magnetically de- 
localized but it is electrically localized. 



(0|aiiA(r-i,0)-D,(r2,0). 
= (0|A(r-i,0)L)j(r2,0), 



,A(r„,0)4|0) 
.A(r-„,0)|0). 



(38) 



The energy operator H of the electromagnetic field, 
unlike the photon number operator [ (|20p ]. is given by a 
local expression 



H= d 



D{r,t)-D{r,t) , B{r,t)-B{r,t) 



2e 



2fi 



(39) 



= ^hc Jd^r[d^{r,t)-d{r,t) + b^{r,t)-b{r,t)] . (40) 

The two terms in the last integrand are not equal to 
the electric and magnetic energy densities that appear 
in Eq. because they do not contain the nondiagonal 
terms that change the total number of photons. These 
terms canceled out only due to the integration over all 
space. However, in a restricted sense the operators 



Heir^t) = dt(r,t)-d(r,t), 
7i:™(r,i) = bt(r,t)-6(r,i). 



(41a) 
(41b) 



can be interpreted as the photon electric and magnetic 
energy densities because the nondiagonal terms do not 
contribute to the expectation values in n-photon states. 
Note that the total magnetic and electric energies are 
equal. 

The expectation values of Tig and 7i„i at t = in a 
one-photon state can be calculated with the use of com- 
mutation relations (|17p . Due to the complementarity be- 
tween the electric localization and magnetic localization, 
the full energy density Tie + "Hm in a one-photon state 
cannot be sharply localized. One can only localize sepa- 
rately either the electric or the magnetic part. 
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The same properties hold for the photon number den- 
sity operator (PT|) . In order to obtain the magnetic local- 
ization we choose the photon state created by operator 
with the wave function <p in the form 



(42) 



The expectation value of the magnetic part of the photon 
number density is then localized in R, 

{0\a4>\r,t)-4>{r,t)a^\0)^\iPn{r)\^ (43) 

However, owing again to the presence of the helicity op- 
erator, the electric part is delocalized, 

{0\aiPHr,t)4{r,t)a^0) = |x^fl(r)|2. (44) 

Analogous results are obtained, mutatis mutandis^ for the 
electric part of the photon number density that can be 
localized while the magnetic part is delocalized. 



VII. DESTRUCTION OF PHOTON 
LOCALIZATION BY TIME EVOLUTION 

Up to this point we have not used the Maxwell equa- 
tions or any other elements of the photon dynamics. All 
conclusions were based on the instantaneous properties 
of the system. The impossibility of a simultaneous elec- 
tric and magnetic localizations follows directly from the 
complementarity between electric and magnetic nature 
of a photon. 

Even this partial (electric or magnetic) localization can 
be achieved only momentarily. It is immediately de- 
stroyed by the time evolution. The solution of the initial 
value problem for the Maxwell equations (see, for ex- 
ample, [l3]) shows that the electric and magnetic field 
operators at t ^ depend on both field vectors at f = 0, 

D{r,t)^ Jd^r'[dtD{r-r',t)D{r',0) 

+ D{r-r',t)V xB{r',0)/fi], (45a) 

B{r,t)^ Jd\'[dtD{r~r',t)B{r',0) 

- D{r-r',t)V X D{r',0)/e], (45b) 
where D{r,t) is the Jordan-Pauli function p^ . 



D{r,t) 



4:TTcr 



[S{ct - r) - S{ct + r)] 



(46) 



Thus, the expectation values of the magnetic fields at t, 

{OlanB, (n , t)Bj (rj, t) . . . Bk (r„ , i)aj^ |0) , (47) 

are coupled to the expectation values involving both the 
magnetic and the electric fields at t = 0. Since a magnet- 
ically localized photon is not electrically localized, even 
its partial localization is an ephemeral effect. 



VIII. DIFFUSE LOCALIZATION OF PHOTON 
STATES 

Since we have shown that sharp localization in a re- 
stricted region is impossible, we can only try to localize 
photons approximately. It turns out that we can con- 
struct photon states m| for which the energy density 
has an exponential tail with an arbitrarily small scale 
factor /. The simplest example of such a state is charac- 
terized by the square-root exponential fall-off. Stronger 
localization is also possible but then we have to resort to 
special functions. Let us choose the wave function in the 
form [H 



tpoi'T) — mlm —m sinK_, (48) 



where m is a constant vector that includes also the nor- 
malization factor and 



K± 



+ (r//)2±l. 



(49) 



With the help of the projector [([S])] we can construct the 
eigenfunction t/'ff (*") = P+'4^d{i') of the helicity operator. 
In Appendix [Cl we evaluate jol'D and obtain a function 
with the same exponential fall-off as ipn, 



r X m e 







/ 




X 


K+ 


cos K_ 






2r + 


/ 2r 



(50) 



The creation operator constructed with the use of the 
function xpn will have the same commutation relations 
(up to factors ±z) with both electric and magnetic fields. 
In this case the electric and magnetic localizations will 
be governed by the same law: the expectation values of 
all the field vectors, electric and magnetic, will fall-off 

when one moves away from the localization 



as e 

center. Therefore, we can choose the scale factor I so that 
at a given distance all field correlation functions will be 
arbitrarily small. Of course, these functions will spread 
out with the speed of light during time evolution [24| . 
Diffuse localization with an exponential tail also holds 
for the photon number density. 



IX. SHARP LOCALIZATION OF COHERENT 
STATES 

All the limitations of the localization apply to the 
states with a given number of photons. Coherent states, 
as we shall show now, are not subject to such limita- 
tions. The easiest way to see this employs the displace- 
ment operator The displacement operator is usu- 
ally constructed from the creation and annihilation op- 
erators that result from a mode decomposition. Here 
we shall construct the displacement operator from our 
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space-dependent operators. The most general coherent 
state can be written in two equivalent forms 

exp Jd^r [ip{r)-d'i{r)- ip*{r)-d{r)] |0) 

= exp J d^r [ifi{r)-b^r) - ip*{r)-b{r)] |0). (51) 

This formula expressed in terms ol the field operators is 

\a,l3)=expifd^r\cx{r)-D{r)+l3{r)-B{r)\ |0). (52) 



The two real vector functions cy.{r) and /3(r) are certain 
combinations of the original complex functions (^(t) or 
ip{r), but at this point we treat them as the primary 
objects. In order to obtain a sharply localized coherent 
state we assume that a. = ckR and j3 = (3r vanish outside 
of some region R. The expectation values of any number 
of field operators in such a localized state are 

(afl, /3fl I a (n , O)^, (r2 , 0) . . . Ck (r„ , 0) | , /3i?) , (53) 

where the operators C stand for D or 13. When all points 
Ti lie outside of i?, we can freely move the operator ap- 
pearing in Eq. (|52p to the left since the field operators 
commute for space-like separations. Next, we use the uni- 
tarity of displacement operator to convert Eq. (|5F 
to the vacuum expectation value, 



(0|Q(ri,0)C',(r2,0)...Cfe(r„,0)|0). 



(54) 



Thus, no measurements in the sense of [ill, |l2| of the 
electric and magnetic fields carried out outside of R will 
reveal the presence of the coherent state localized in R. 
Coherent states may be sharply localized in an arbitrarily 
small region. The localization region remains sharp but 
it will expand in time with the speed of light, as seen in 
formula (|i5|) . 

The localization of the field operators implies the lo- 
calization of the energy density. One can also obtain the 
full localization of the photon number density by choos- 
ing the functions Q;(r) and (3{r) in form (|^^ . 



at an instant. This localization is fragile, it is immedi- 
ately destroyed by time evolution. Due to lack of a sharp 
photon localization, we constructed diffuse photon states 
that are simultaneously electrically and magnetically lo- 
calized with finite but arbitrary accuracy. We have also 
shown that there are no restrictions on the localization 
of coherent states. 
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APPENDIX A 



In our calculations we need the following special cases: 



dh 



I 



d^r 



167r 



1^' _ j.|5/2|^ _ ^//|5/2 |^/_^"|2' 



dh 



An' 



(Al) 



(A2) 



(A3) 



_ .^|5/2|^ _ ^//|3/2 _ .^//| ' 

of the general formula 

■Ka , ,7r/? . a + P ^ 



\r' — r\"\r — r"\P 



TTm = sin(— )sin(— )sin7r(— — ) 



8T{a + 13- 4)r(2 - a)r{2 - (}) 

1^/ _ .^//|ct+/3-3 



(A4) 



X. CONCLUSIONS 



APPENDIX B 



We have shown that the symmetry between electricity 
and magnetism permeating the classical and quantum 
theories of the electromagnetic field enables us to split 
the notion of photon localizability into two separate no- 
tions: the electric and magnetic localizabilities. We de- 
fined operationally the electric (or magnetic) localization 
of photon states as follows: the electric (or magnetic) 
footprints of electrically (or magnetically) localized pho- 
tons are not detectable outside of the localization region. 
Since the electric and magnetic fields are complemen- 
tary, electric and magnetic localizabilities are mutually 
exclusive. We have further shown that even this par- 
tial sharp localization of a photon can be achieved only 



We give here a general proof by reductio ad absurdum 
that a simultaneous sharp localization of both ^^}{r) and 
Xi/'(t) is not possible. 

Take a function xljR{r) sharply localized in R. Let 
us suppose that the function X*}}R{r) is also sharply lo- 
calized in some region R' . Then the sum of these two 
functions ip+{r) = (1 + x)'^B.{f) is a function with com- 
pact support in the region RU R' . Due to the pres- 
ence of the factor {1 + x), it is also an eigenfunction of 
the helicity operator belonging to the eigenvalue 1. The 
Fourier transform of every integrable function with com- 
pact support is an entire function (see, for example, [2^). 
The Fourier transform ip+{k) of the function ip+{r) is an 
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eigenfunction of the hehcity operator. Therefore, accord- 
ing to Eq. ([5]), it must obey the relation 



k X ip+{k) 



1,2 



1,2 



kl ^+{k). 



(Bl) 



In the first step we used the formula (Eq. 3.472.5. in [26| ) 



'^^e-f--'/- = (-1)" /^|le-2v^. (C3) 



We arrived at a contradiction because the left hand side is 
an entire function but the right hand side is not because 
it contains a square root. Therefore, the function xxpf{{r) 
cannot be sharply localized. 

A less formal but more physical argument goes as fol- 
lows. According to the formula ([7]), the function xipnir) 
has a l/r^ tail extending to infinity. This fall-off is 
generic but we could obtain a faster decay by fine-tuning 
the function xpii{r). Namely, we can require the van- 
ishing of its successive moments. However, in this way 
we can never make the function xipR{r) to vanish com- 
pletely outside of R because ipRir) cannot have all the 
moments equal to zero. 



APPENDIX C 

We calculate the function xy_D using the expression 
for the convolution in terms of the Fourier transforms. 



d'r'fir-r')gir') 



(2^ 



e^'^-^mgik). (CI) 



The Fourier representation of the function appearing in 
Eq. ([HI) has the form 



Ime-Vi-^-^/' _ 1 /-^ sin(fcr) 



-kl-l/kl 



ryfid Jo '^^ k^'^ 
27r3/2 f d'^k „-^^e-'='-i/fci 



VI J (2^)3 fc5/2 



(C2) 



The imaginary parts in Eqs. (jC2p and (jC6|) can be ex- 
tracted with the use of the following algebraic identity: 



V2Va + ib^ ^/\a + ib\ + a + J^VW + ib\ - a . (C4) 



Putting together formulas (|Cip . and (|C2|) . we ob- 
tain 



XVoir) = V X Jd^r 



2TT^\r — r' 



27r3/2 



Vx /^e^^-- 



-kl-l/kl 



(27r)3 



1 ^ 1 Z"^^, sin(fcr) _ki-i/ki 
V X - / dk^^jj^e ^/'''m 







■kI 

1 I f°° Q-k(l-ir)-l/kl 

V X Im- / dk -^7- m. 

Using again Eq. (jC3[) . we arrive at the final result 



(C5) 



r X m. 



Tm 



2r r V I 



(C6) 
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